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Abstract. We present STAR's measurements of directed flow for charged hadrons 
in Au+Au and Cu+Cu collisions at y'SNN = 200 GeV and 62.4 GeV, as a function 
of pseudorapidity, transverse momentum and centrality. We find that directed flow 
depends on the incident energy, but not on the system size. We extend the validity of 
limiting fragmentation hypothesis to different collision systems. 



1. Motivation 

Directed flow is quantified by the first harmonic (v{) in the Fourier expansion of the 
azimuthal distribution of produced particles with respect to the reaction plane [lj. It 
describes the collective sideward motion of produced particles and nuclear fragments, 
and carries information from the very early stage of the collision [2]. Charged-hadron 
Vi results at RHIC have been previously reported over a wide range of incident 
energies in Au+Au collisions [31 HI El E], and they support the limiting fragmentation 
hypothesis [U El II]- In RHIC run V (2005), lighter nuclei were collided at a/snn = 200 
GeV and 62.4 GeV, which enables us to study the system-size dependence of directed 
flow, as well as the incident-energy dependence. 



2. Data sets and the study 

This study is based on eight million 200 GeV Au+Au events, five million 62.4 GeV 
Au+Au events, twelve million 200 GeV Cu+Cu events, and eight million 62.4 GeV 
Cu+Cu events. All of them were obtained from a minimum-bias trigger. All errors are 
statistical. Tracks of charged particles are reconstructed by STAR's main TPC [8] and 
forward TPCs [9], with the pseudorapidity (77) coverage of \r]\ < 1.3 and 2.5 < \r)\ < 4.0, 
respectively. The centrality definition and cuts used are the same as in Ref. [TO] . 

At RHIC energies, it is a challenge to measure V\ accurately due to the small 
signal of V\ itself and the large systematic error arising from non-flow correlations. In 
this talk, STAR's preliminary v i{ZDC-SMD} results are obtained with non-flow effects 

I For the full list of STAR authors and acknowledgements, see appendix 'Collaborations' in this volume. 
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Table 1. The resolution of the lst-order full event plane provided by STAR ZDC- 
SMDs, as determined from the sub-event correlation between east and west SMDs. 
The errors in the table are statistical. 



Centrality 


200 GeV Au+Au 


62.4 GeV Au+Au 


200 GeV Cu+Cu 


62.4 GeV Cu+Cu 


70% - 80% 


0.296 ±0.003 


0.179 ±0.005 






60% - 70% 


0.348 ±0.003 


0.185 ±0.004 






50% - 60% 


0.382 ±0.002 


0.176 ±0.005 


0.140 ±0.003 


0.039 ±0.011 


40% - 50% 


0.397 ±0.002 


0.167 ±0.005 


0.144 ±0.003 


0.043 + 0.011 


30% - 40% 


0.390 ±0.002 


0.138 ±0.006 


0.147 ±0.003 


0.032 ±0.015 


20% - 30% 


0.365 ±0.002 


0.110 ±0.008 


0.121 ±0.004 




10% - 20% 


0.309 ±0.003 


0.081 ±0.010 


0.095 ±0.005 




5% - 10% 


0.220 ±0.006 




0.059 ±0.008 




0-5% 


0.120 ±0.002 




0.059 ±0.008 





minimized by determining the event plane from spectator neutrons, the same approach 
as used in Ref. [5], [6]. The resolution, as defined in Ref. pQ, of the first-order event 
plane reconstructed with STAR ZDC-SMDs [TT] is listed in Table [TJ The magnitude of 
the event plane resolution increases with the spectator v i and the detector efficiency of 
the ZDC-SMDs, and the latter is smaller for a lower incident energy and/or a smaller 
collision system. The method fijZDC-SMD} fails in some centralities, where the event 
plane resolution is consistent with zero. 

3. Results 

Figures in this talk follow the same convention as in Ref. [5]. 

Fig. [1] shows v\ as a function of pseudorapidity, r], in 30% — 60% most central 
Au+Au and Cu+Cu collisions at a/snn = 200 GeV and 64.2 GeV. Data points in Fig. [T] 
fall into two bands, and the higher band corresponds to the lower energy, which reveals 
the incident-energy dependence of viiji). For the pseudorapidity and centrality range 
we studied, vi(r]) in Au+Au collisions and Cu+Cu collisions are consistent within errors 
if both of them were obtained at the same incident energy. It is interesting that v i does 
not change from Au+Au collisisons to Cu+Cu collisions, while the system size is reduced 
by 1/3. A good consistency of v± between Au+Au and Cu+Cu collisions is observed 
even for the region near midrapidity, where v<i in Cu+Cu collisions is considerably lower 
than that in Au+Au collisions [12] . Unlike v^/e which scales with \^^ L (interpreted to 
be the participant density [12] or the system length [13]), V\ is found to be independent 
of the system size. Instead, it scales with the incident energy. A possible explanation 
to the different scalings for t> 2 /e and V\ might come from the way in which they are 
developed : to produce t>2, intensive momentum exchanges among particles are needed 
(and the number of momentum exchanges is related to the participant density or the 
system length), while to produce one in principle needs only different rapidity losses 
(related to the incident energy) for particles with different distances to the center of the 
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Figure 1. Charged-hadron vi vs. r/, in Au+Au Figure 2. Charged-hadron i>i vs. ?y — j/bc 



and Cu+Cu collisions at 200 GeV and 62.4 GeV. 
The plotted errors are statistical. 



Au+Au and Cu+Cu collisions at 200 GeV and 62.4 
GeV. The plotted errors are statistical. 



participant zone. 

In Fig. |2l we re-arrange the data points of Fig. [I] in the projectile frame, where zero 
on the horizontal axis corresponds to the beam rapidity, ybeam, for each of the incident 
energies. Within three units from ?/beam, most data points fall into a universal curve 
of V\ as a function of r\ — ?/bcam- This incident-energy scaling behavior of directed flow 
suggests that the limiting fragmentation hypothesis [7j holds even for different collision 
systems. 

The transverse-momentum (p t ) dependence of v\ for charged hadrons is shown in 
Fig. [3j The upper panel shows the results measured in the main TPC, and the lower 
panel, in the forward TPCs. Since Vi(j],pt) is asymmetric about r] — 0, the integral 
of Vx(r),pt) over a symmetric r] range goes to zero. We change vi(r],p t ) of particles 
with negative rj into —Vi(—i],p t ), and integrate over all r]. At 200 GeV, V\{p t ) in 
Au+Au and Cu+Cu collisions are very close to each other (better visible in the forward 
pseudorapidity region), and V\{jp t ) saturates above p t ~ 0.8GeV/c. At 62.4 GeV, V\{jp t ) 
reaches its maximum in magnitude at p t w 0.6 GeV/c, then saturates (\r]\ < 1.3) 
or slightly decreases (2.5 < \n\ < 4.0). Within statistical errors, the incident-energy 
scaling holds in Vi(p t ) as well. 

Fig. H] shows charged-hadron ^-integrated v\ as a function of centrality for the 
two collision systems at 200 GeV and 62.4 GeV, and the results are presented in two 
panels defined in the same way as Fig. [3] On the whole, the magnitude of integrated 
V\ approaches zero in central collisions and increases as the collisions become more 
peripheral. At 200 GeV, the pt-integrated v\ in Au+Au collisions is consistent with that 
in Cu+Cu collisions in almost all centrality bins, in regions both near the midrapidity 
(|?7| < 1.3) and away from it (2.5 < \r]\ < 4.0). This system-size independence is 
seemingly true at 62.4 GeV. We can see that the incident-energy scaling of v\ works 
well as a funtion of centrality, instead of the participant density, with which the i^/e 
scales [12]. This indicates that the initial geometrical shape, but not the density of 
produced particles, plays an important role in the formation of directed flow. 
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Figure 3. Charged-hadron v\ vs. p t , in 30%— 60% 
Au+Au and Cu+Cu collisions at 200 GeV and 62.4 
GeV. The upper panel shows the results measured 
in the main TPC, and the lower panel, in the 
forward TPCs. The plotted errors are statistical. 



Figure 4. Charged-hadron p t -integrated V\ vs. 
centrality, in Au+Au and Cu+Cu collisions at 200 
GeV and 62.4 GeV. The upper panel shows the 
results measured in the main TPC, and the lower 
panel, in the forward TPCs. The plotted errors are 
statistical. 
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